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[~^. ' We present a consistent regularization procedure for calculating hadron structure func- 

tions in a bosonized chiral quark model. We find that the Pauli-Villars regularization 
scheme is most suitable. We also summarize the phenomenology of structure functions 
calculated in the valence quark approximation. 
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^^ ■ 1 Introduction 

in 

o 



The strong interaction contributions to the cross-sections in deep inelastic scat- 
tering are parameterized by the hadronic tensor W^ v . The associated form factors 
^^ ■ depend on the Lorentz-scalars Q 2 = —q 2 and x = Q 2 /2p ■ q, with p and q being 

""p^ | the momenta of the nucleon and the exchanged photon, respectively. The most 

convincing evidence for the quark-substructure of the nucleon stems from the so- 
called scaling behavior of these form factors implying that certain combinations 
j~< thereof i.e. the structure functions, depend on Q 2 only logarithmically. Without 

bound-state wave functions in continuum QCD it has proven difficult to calcu- 
late the hadron structure functions from first principles. Fortunately the operator 
product expansion supplies an avenue to calculate nucleon structure functions: To 
leading twist (expansion in l/Q 2 ) the moments of the structure functions factor- 
ize into non-perturbative and perturbative pieces. The former are tractable within 
effective low energy theories for QCD while the latter represent the logarithmic 
scaling violations which can be computed utilizing the renormalization group equa- 
tions M. Alternatively one may first sum the series i.e. calculate the leading twist 
structure functions at the "model scale", Qq, and subsequently incorporate the 
radiative correction using the DGLAP evolution procedure. 

In a field theory calculation the most suitable starting point to compute hadron 
structure functions is the forward Compton amplitude for hadron photon scattering 



T^(p,q;8)=iJd?Se**(p t s\T{Jr(t)4(p)}\p,8), (1) 
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since the time-ordered correlation function of electro-magnetic currents, J M can 
be calculated in the functional integral formalism. In turn the hadronic tensor is 
obtained from its absorptive contribution, 



1 
2^ 



Wm„ (p, q; s) = — Im T M „ (p, g; s) . (2) 



The structure functions are then extracted from W 7 ^ by adopting the Bjorken limit, 
Q 2 — > oo with a; fixed. 



2 Generalities on the Structure Functions in the NJL— Model 

The results and discussions presented in this section are based on ref. 0. The 
interaction in the Nambu-Jona-Lasinio (NJL) model is a chirally invariant four 
fermion coupling, with coupling constant G. Employing functional integration tech- 
niques we may recast it into an effective meson theory |3]|. What remains is to 
properly regularize the UV divergence of the resulting determinant of the Dirac 
operator in the background of the meson fields. We utilize the Pauli-Villars (PV) 
regularization scheme which provides a formulation of the bosonized NJL model 
completely in Minkowski space [Q . This is suitable when applying Cutkosky's rules 
to extract the hadronic tensor via (g) . In addition the P V scheme yields the proper 
scaling of the structure functions since the UV cut-off is additive in both the loop 
momenta and the (constituent) quark masses. 

In the NJL chiral soliton model regularizing the nucleon structure functions 
has only recently been addressed in a satisfactory manner []2| . In the first calcula- 
tions |j |(j, 0i contributions from the polarized vacuum to the structure functions 
were approximated within the (constituent) valence quark approximation. This has 
been motivated by the observation that the regularized vacuum contributions to 
static nucleon properties are small once the self-consistent soliton solution is deter- 
mined. Later approaches including vacuum contributions imposed a single subtrac- 
tion regularization scheme at the level of (constituent) quark distributions M, M. 
The regularization attributed to a given structure function was conjectured to equal 
that of the corresponding sum rule associated with static nucleon properties. This 
approach is obviously inadequate for structure functions whose sum rules are not 
associated with a matrix element of a local operator. Also, the single subtraction 
scheme leaves the gap-equation determining the VEV, m = (S) undefined. 

To calculate the electro-magnetic current operator appearing in the Compton 
tensor ([[]) we introduce sources v^ which are conjugate to this current 

iD = i# - (S + ijnP) + Qj =: iD (,r) + Qi> (3) 

iD 5 = -i# -(S- if S P) - Qi =: iD^ 3 - Qj, , (4) 

with Q = diag(2/3, 1/3) denoting the quark charge matrix. S and P are respectively 
scalar and pseudo-scalar matrix valued fields. In addition to the common Dirac 
operator D, is the 75 -odd Dirac operator D5, which in Euclidean space would be 
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associated to the Hermitian conjugate of D. The bosonized and PV regularized 
form of the NJL-model then reads ( ^\ a — ^T CiA 2 = 0, Ao = 0, co = l) , 



-4njlK] = i d 4 xtY[S 2 + P 2 + 2m S] 

2 

i^f V cTrlog [ DD 5 + A 2 - ie] - i^Tr log [-D (Ds)" 1 



(5) 



- 1 '2 ^ 

i=0 



The fact that the last part of (g) is not regularized ensures that the ABJ-anomaly is 
completely recovered. The time ordered matrix element of the currents is obtained 
by expanding the action to quadratic order in the sources v^. In ref. (2| we have 
shown that considerable simplifications arise in the Bjorken limit; in particular the 
propagator between two insertions of the sources becomes that of a free quark, 
yielding the action 

i— n ^ J 



1=0 



!^£tv J LnWnW 



+ ___Tr { ( -D W D 5 



(6) 



The above mentioned decomposition of the action into regularized and unregular- 
ized pieces has enforced the specification {S^ pva = g w g va + .9^.9^ - g^g P <j), 

lyTfplv = S wv <j^ - ie lipva j (T 'Y 5 while (7^7/^)5 = S^^-y' 7 + ie P p Va -y a ^ 5 . (7) 

The latter reduction in (R) which is reflected in that of (^ ($) 1^)5, arises in order 
that the regularization of the structure functions be consistent with the corre- 
sponding sum rules)^)]. Finally from (g) we calculate the bilocal current correlation 
function and apply Cutkosky's rules to obtain the absorptive part (0). 

The pion fields are identified as the angular variables on the chiral circle, 
S + iP = mexp(i<?T • if/m) with the constituent quark mass being the vacuum 
expectation value (VEV) of the scalar field, m = (S). The quark-pion coupling 
constant is extracted from the normalization of the pion field yielding 

T* =ANc lalf dX ^ Ii{q2 ^\i .' (8) 

g «<? Jo q 2 =ml 

where 

n( 9 2 ,x) = -i^Ci^ [-k 2 -x(l-x) q 2 +m 2 +A*-ie]- 2 (9) 

4 =o [Z7T) 

is the pion-polarization function parameterizing the quark loop. Upon expanding 
D' 7r ) and Dg up to quadratic order in the pion fields we obtain the virtual Comp- 
ton amplitude for pion-photon scattering and extract the pion structure function 

1 ) This issue is not unique to the PV regularization scheme; in fact all schemes which regularize 
the sum log (D) + log (D5) but not the difference log (D) — log (D5) require this treatment. 
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to be 



F(x) 



^^[flW,*)] 



(10) 



In the chiral limit, m^ = 0, the pion structure function turns out to be independent 
of x and simply reproduces the pion-quark coupling (|S|) . Although this result has 
been known [|10| it serves to verify the approach defined by the expansion (|6|). 

In this model the nucleon emerges as a chiral soliton pl| . The starting point 
for constructing the static soliton is the Dirac Hamiltonian, h is defined from (J|), 



iD w = 0(id t - h) and iT> 



(t) 



(-id t -h)p. 



(11) 



It is diagonalized for the well-known hedgehog configuration of the background 
meson fields, 



Zi'I'a = e Q $ Q with h = a ■ p + m/3exp I if ■ T'fsQfr] 



(12) 



yielding eigen-spinors *& a and energy eigenvalues e a . For this mesonic configura- 
tion the action functional is expressed as a regularized sum over the eigenvalues e a . 
This yields an energy functional, E[@] which equals, to leading order in 1/Nc, the 
nucleon mass. The minimization of E[Q] determines the localized profile function 
0(r), self-consistently. Nucleon states are generated by introducing collective co- 
ordinates for the position of the soliton as well as its orientation in coordinate and 
flavor spaces. Canonical quantization of these coordinates yields the nucleon states 
which enter the matrix elements (flj). Finally the structure functions are obtained 
by contracting the corresponding hadronic tensor (0) with pertinent projectors. As 
an example we list the resulting expression for the polarized structure function g\ , 



9i{x) 



.M N [e]N c fdtd 



:-!() 



2tt 



?/*/ 



dX 

2^ 



iM N [Q]x\ 



(13) 




UT 



Kj 



[®i(£)T3(l-a 3 )j 5 y a (£+)e 



30 



a 2 [Q 

[*1(£> 3 (1 



4 



N 



■i.LuX 



P 

— ILOX 




a 3 )75*a(^) e 
(uj + e a )(u + ep) 



N 



^(l)T 3 (l-a 3 )75* Q (^)e^ A 



<(0r 3 (l-«3)75*o(C)e 11 



ie)(cu 2 



•te) 



ar 3 



*]»(!) (1 -oa) 75* Q (el)e~^ A + * J(I5 (1 -« 3 ) 75* Q (C)e^ A 
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where £ A = £ ± A<§3 label shifted coordinates. The subscript 'p' indicates that 
the spectral integral is restricted to those values of u> which cause the respective 
denominators of the spectral functions to vanish. Similar expressions are obtained 
for /i and <?2 [pi ■ For the polarized structure functions is important to remark that 
the isoscalar contribution, which is next to leading order in 1/Nq does not undergo 
regularization while the leading order isovector contribution does. As in the study of 
static nucleon properties the NJL model we find that either the isoscalar or isovector 
contributions to any structure function become regularized, but not both. The result 
that the isovector part is regularized was anticipated because we know how to relate 
the zeroth moments to nucleon charges |2| . With regard to the unpolarized structure 
function fi{x) we find the reverse situation: a regularized leading order isoscalar 
and an unregularized next to leading order isovector contribution. The latter piece 
enters the Gottfried sum rule implying that it is unregularizcaQ- 

The verification of the sum rules, which relate moments of the structure func- 
tions to static properties of hadrons, is lengthy but straightforward. Essentially it 
can be traced back to the expansion m) which can be formally interpreted as a 
first order expansion in the operator Q 2 (^(^) 1^)5. Such first order expansions 
are related to static properties. 

3 Numerical Results for the Nucleon Structure Functions 

Numerical results are still unavailable for the structure functions e.g. (O) as 
projected from the consistently regularized hadronic tensor. Thus, we confine our- 
selves to presenting results from the valence quark approximation H, pi . To compare 
the model structure functions with data, we calculate at the intrinsic scale of the 
model in the nucleon's rest frame (RF), transform them to the infinite momen- 
tum frame (IMF) and evolve the results to scales commensurate with the data. 
We note that the Bjorken limit singles out the IMF which is preferred for structure 
functions calculated from extended field configurations such as the chiral soliton be- 
cause in this frame this class of structure functions possesses proper support [ |12| . 
The DGLAP evolution determines this scale parameter Qq at which the model 
approximates QCD. 

The unpolarized structure functions are calculated from the symmetric portion 
of the hadronic tensor, W^ v + W V)1 . In Fig. 1 we display the result of the structure 
functions entering the Gottfried sum rule for e—N scattering as calculated in the RF 
and in the IMF. The DGLAP evolution determines our scale to be, Q\ « 0.4GeV 2 . 
The resulting structure function reproduces the gross features of the experimental 
data. In addition depending on the other model parameter to, the Gottfried sum 
rule, Sg = J dx (F^ — F£ n ) jx ranges between 0.26 — 0.29 which exhibits the 
desired deviation from the naive quark model prediction of 1/ 3 and is in reasonable 
agreement with the experimental range, 0.235 ± 0.026 p3[. 



2 ) This is in contrast to previous studies fej, B| where this sum rule was regularized based upon 
the assumption that it behaved equivalently to the Adler sum under regularization. 
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Fig. 1. Left Panel: The unpolarized structure function entering the Gottfried sum rule: 
RF, IMF and leading order QCD evolution to Q 2 = 3GeV 2 , with full 1/JV C contribution; 

see M. Right Panel: gi(x). 




5.0GeV . Right Panel: Chiral odd structure function \iP T [x,Q 



The spin-polarized structure functions g± and g 2 are extracted from the anti- 
symmetric combination, W^ — W vil . Results are shown in Fig. 2 and are compared 
with the data from SLAC |14J. In particular we note the sizeable twist two and 
three contributions. The twist two portion is given by the Wandzura-Wilczek con- 
tribution, the first two terms in 



52 (x, Q 2 ) = - 9l (x, Q 2 ) + I -^-g x (x, Q 2 ) + g 2 (x, Q 2 ) 

J x y 



(14) 



where g 2 denotes the twist three contribution. In the evolution of the twist- 2 portion 
we have restricted ourselves to the leading order in a s since the next to leading 
order evolution of the twist-3 ~g 2 \ x iQ) i s n °t known (The twist-3 piece is known 
only in the large 1/N C limit fig).). In ref. H the Bjorken J dx{g\{x)—gi{x)) = gA/6 
and Burkhardt-Cottingham J dx g 2 (x) = sum rules have been explicitly verified. 
Having calculated the spin (un)polarized structure functions, /i and g\ and 
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their constituent quark components we investigate the chiral odd flavor distribu- 
tions and their charge weighted average nucleon "structure functions" |l6[ . Here we 
confine ourselves to the presentation of results relevant in the context of the Soffer 
inequality J17J. This inequality relates the nucleon chiral odd quark distribution 



functions to both the unpolarized /{ 
counter-parts. 

» 



O) 



(x, Qq) and polarized g\ (x, Qq) chiral even 



ti u> (x,Q 2 )+g[ u >(x,Q 2 )>2h^(x,Q 2 ) 



(15) 



The superscript refers to the flavor combination which projects onto up-quark 
quantum numbers. Note again, that this projection refers to the constituent quarks 
which contain some non-perturbative gluonic distributions. Fig. 3 clearly demon- 
strates that the inequality is satisfied at the model scale. 





Fig. 3. Left Panel: The Soffer inequality for the chiral even combination in the IMF (solid 

line) for the up-quark distributions and the chiral odd structure function 2/i^ {x,Qo) 

(long-dashed line) . Right Panel: Down quark distributions. 



4 Conclusions 

We have summarized our results on the consistent regularization scheme in cal- 
culating hadron structure functions in the bosonized NJL chiral soliton model where 
baryons emerge as chiral solitons. We have calculated the Compton amplitude in 
the Bjorken limit from the PV regularized action functional with the constraint of 
preserving the anomaly structure of QCD. Introducing pertinent external sources 
enables the unambiguous identification of the quark currents which in turn allowed 
us to calculate the bilocal current correlation function entering the forward Comp- 
ton amplitude. This scheme satisfies the leading order scaling properties of nucleon 
structure functions. The sum rules for both polarized and unpolarized are fulfilled 
because in the associated moments Cutkosky's rules extract - similar to the Cauchy 
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integrals - the matrix elements associated with static nucleon properties H . These 
results represent the first step towards a fully consistent calculation of structure 
functions in the NJL model. We have also summarized our results in the valence 
quark approximation. Due the sizeable contribution from this state to all nucleon 
properties in the NJL model, not surprisingly, this approximation reproduces the 
general trends of the data on structure functions. 
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